A Monte Carlo Synthetic Acceleration Method for the Non-linear,time-dependent Diffusion Equation
نویسنده
چکیده
We present a Monte Carlo synthetic-acceleration method for solving the time-dependent, nonlinear, equilibrium diffusion equation in three-dimensions. The new scheme uses the adjoint Monte Carlo method as a relaxation step that accelerates standard Jacobi iteration. Results show that this method is 40% faster than regular Jacobi iteration and is competitive with Jacobi-preconditioned Conjugate Gradient methods. Furthermore, the new method is not limited to symmetric, positive-definite systems, and therefore, it can be used for non-symmetric systems. Such systems arise in fully coupled, nonlinear-consistent (Newton) solvers.
منابع مشابه
A Monte Carlo simulation technique for assessment of earthquake-induced displacement of slopes
The dynamic response of slopes against earthquake is commonly characterized by the earthquake-induced displacement of slope (EIDS). The EIDS value is a function of several variables such as the material properties, slope geometry, and earthquake acceleration. This work is aimed at the prediction of EIDS using the Monte Carlo simulation method (MCSM). Hence, the parameters height, unit specific ...
متن کاملPlanar and SPECT Monte Carlo acceleration using a variance reduction technique in I131 imaging
Background: Various variance reduction techniques such as forced detection (FD) have been implemented in Monte Carlo (MC) simulation of nuclear medicine in an effort to decrease the simulation time while keeping accuracy. However most of these techniques still result in very long MC simulation times for being implemented into routine use. Materials and Methods: Convolution-based force...
متن کاملGoodness of fit of stochastic differential equations
We propose a method to test for lack-of-fit of an estimated stochastic differential equation. The method is based on Monte Carlo simulation of trajectories between neighbour observations and, thus, it does not rely on the availability of explicit expressions of the conditional densities. Consequently, both non-linear models and models with state-dependent drift and diffusion can be handled. The...
متن کاملThe Use of Monte-Carlo Simulations in Seismic Hazard Analysis in Tehran and Surrounding Areas
Probabilistic seismic hazard analysis is a technique for estimating the annual rate of exceedance of a specified ground motion at a site due to the known and suspected earthquake sources. A Monte-Carlo approach is utilized to estimate the seismic hazard at a site. This method uses numerous resampling of an earthquake catalog to construct synthetic catalogs to evaluate the ground motion hazard a...
متن کاملDiffusive Acceleration of Ions at Interplanetary Shocks
Heliospheric shocks are excellent systems for testing theories of particle acceleration in their environs. These generally fall into two classes: (1) interplanetary shocks that are linear in their ion acceleration characteristics, with the non-thermal ions serving as test particles, and (2) non-linear systems such as the Earth’s bow shock and the solar wind termination shock, where the accelera...
متن کامل